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Recent work on ultracold polar molecules, governed by a generalization of the t-J Hamiltonian, 
suggests that molecules may be better suited than atoms for studying d-wave superfluidity due to 
stronger interactions and larger tunability of the system. We compute the phase diagram for polar 
molecules in a checkerboard lattice consisting of weakly coupled square plaquettes. In the simplest 
experimentally realizable case where there is only tunneling and an XX-type spin-spin interaction, 
we identify the parameter regime where d-wave superfluidity occurs. We also find that the inclusion 
of a density-density interaction destroys the superfluid phase and that the inclusion of a spin-density 
or an Ising-type spin-spin interaction can enhance the superfluid phase. We also propose schemes 
for experimentally realizing the perturbative calculations exhibiting enhanced d-wave superfluidity. 

PACS numbers: 67.85.-d, 71.10.Fd, 74.20. Mn, 74.20.Rp 



I. INTRODUCTION 

The Hubbard Hamiltonian is believed to contain some 
of the ingredients necessary to explain high-temperature 
superconductivity in cuprates [1-3]. The difficulties of 
analytically understanding the Hubbard Hamiltonian in 
more than one dimension suggest the use of experimen- 
tal quantum simulators to investigate the physics of this 
model, and important experimental progress in this di- 
rection has been made with lattices of ultracold atoms 
[4-10] . In the context of high-temperature superconduc- 
tivity, the most relevant regime of the Hubbard Hamilto- 
nian is the limit of strong on-site interactions, in which 
the model reduces to the t-J Hamiltonian [1-3]. Here t 
is the nearest-neighbor tunneling amplitude and J is the 
nearest-neighbor exchange interaction originating from 
second-order virtual hopping. Unfortunately, the ex- 
change interactions are so small that it is extremely dif- 
ficult to observe the associated physics in the cold atom 
implementation [8]. 

Recently, the polar molecules KRb and LiCs have been 
cooled to their electronic, rotational, and vibrational 
ground states [11-14], and KRb has been loaded into 
a three-dimensional optical lattice [15]. This system can 
be used to implement lattice Hamiltonians based on ro- 
tational states of polar molecules [16-32]. Specifically, 
two rotational states of the molecules can be used as 
an effective spin- 1/2 degree of freedom, while dipole- 
dipole interactions mediate "spin" -dependent coupling 
between molecules. For molecules on neighboring sites, 
the strength of these dipole-dipole interactions is ~ 1kHz 
for KRb and ^100 kHz for LiCs. These interactions are 
much stronger than the exchange interactions between ul- 
tracold atoms (<C 1 kHz [9]). Therefore, polar molecules 
seem to be better candidates for the simulation of certain 
condensed matter phenomena. 

Recently, it was shown that polar molecules in optical 
lattices can be used to simulate a highly tunable general- 
ization of the t-J Hamiltonian, referred to as the t-J-V- 



W Hamiltonian [30, 31]. The latter differs from the t-J 
Hamiltonian in the following aspects: it has anisotropic 
XXZ spin-spin interactions Jj_ and J z , an independent 
density-density interaction V , a spin-density interaction 
W, and the interactions are long-range dipolar rather 
than nearest-neighbor. This Hamiltonian is highly tun- 
able, and the strengths of these interactions can, in prin- 
ciple, be varied independently in experiments. In par- 
ticular, the regime J > t can be achieved, which is not 
possible with cold atoms where J originates from second- 
order virtual hopping. Furthermore, in Ref. [31], it was 
shown that, in one-dimension, the t-J-V-W model can 
support enhanced superfluidity relative to the standard 
t-J model. 

In this paper, we use the tunability of the t-J-V-W 
Hamiltonian to find parameter regimes supporting robust 
d-wave superfluidity in one- and two-dimensional systems 
of weakly coupled plaquettes [33-38]. Furthermore, we 
demonstrate that this solvable limit and the associated d- 
wave superfluidity are experimentally realizable. Finally, 
we believe that this limit can provide qualitative guidance 
for the case of the homogeneous two-dimensional lattice. 

Throughout the paper, we consider, for simplicity, an 
average filling of three molecules per plaquette. In the 
simplest experimentally realizable case where J z = V = 
W = 0, referred to as the t-J± Hamiltonian, we find three 
phases: a d-wave superfluid of bound holes, a checker- 
board solid of alternating plaquettes of bound holes, and 
phase separation of bound holes. We find that the ad- 
dition of a density-density interaction proportional to V 
destroys the superfluid phase in the perturbative limit. 
We also find that an Ising spin-spin interaction propor- 
tional to J z or a spin-density interaction proportional to 
W can enhance the superfluid phase for certain parame- 
ter regimes. 

The remainder of the paper is organized as follows. 
In Sec. II, we introduce the t-J-V-W Hamiltonian and 
briefly discuss its experimental implementation with po- 
lar molecules in optical lattices. In Sec. Ill, we ana- 
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FIG. 1. (Color online) Geometry of plaquettes. J schemat- 
ically denotes all dipolar interactions J±, J z , V, and W. 
Within each plaquette, solid blue lines contain tunneling of 
amplitude t and dipolar interactions of strength J, while 
solid red lines only contain dipolar interactions of strength 
J/(2\/2)- Nearest-neighbor plaquettes are then linked with 
tunneling of amplitude t' and dipolar interactions of strength 
J', shown as dotted blue lines, where t' <C t and J' <C 
J. Furthermore, nearest-neighbor and next-nearest-neighbor 
plaquettes are linked with dipolar interactions of strength 
J'/(2y / 2), shown as dotted red lines. The numbering of sites 
in a single plaquette is shown in black. 



lyze the t-J± Hamiltonian in detail solving the single- 
plaquette Hamiltonian exactly and calculating the phase 
diagram perturbatively. In Sec. IV, we analyze the ef- 
fects of the J z , V, and W terms on the superfluid phase. 
In Sec. V, we present proposals for experimentally re- 
alizing the perturbative calculations for one and two- 
dimensional systems of plaquettes. Finally, in Sec. VI, 
we present conclusions. Appendix A gives a brief sum- 
mary of the group theoretic techniques used to solve the 
single-plaquette Hamiltonians exactly and describes the 
symmetries of the solutions. Appendix B gives explicit 
expressions for the basis vectors of the irreducible rep- 
resentations discussed in Appendix A, and Appendix C 
gives explicit expressions for important ground-state en- 
ergies and states. 
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, (1) 



c\, a is the creation operator for a hardcore 
fermionic molecule at the lattice site r with effective 



spin cr, n ra = cl a c ra , n r = n r ^ + n r ±, S+ = c^c^, 

S~ = {S+)\ and S z = {n rt -n rl )/2. The tunnel- 
ing amplitude t rr > and the dipolar interaction strengths 

are t, J±, J Zl V, and W re- 
spectively if the sites r and r' are in the same plaquette 
and are t', J' ± , J' z , V , and W respectively if the sites are 
in neighboring plaquettes. The () signify that the sums 
are taken over nearest-neighbor bonds and the (()) sig- 
nify that the sums are taken over next-nearest-neighbor 
(diagonal) bonds. The next-nearest-neighbor bonds have 
a factor of 1/(2V2) since they are a factor of v2 longer 
than the nearest-neighbor bonds and since the dipolc- 
dipole interaction strength falls off inversely as distance 
cubed. Eq. (1) omits the energies of the states |f) and \l) 
since we work at fixed numbers of up and down molecules. 

The Hamiltonian given by Eq. (1) could be experimen- 
tally realized by loading ultracold polar molecules into an 
optical lattice and applying an external dc electric field 
perpendicular to the plane of the lattice [30, 31]. Two 
rotational states |mo) and \mi) of a molecule form the 
effective spin states |f) and \\), respectively. The prepa- 
ration of these states is discussed below in Sec. VA2. 
Due to the dc electric field, these states have permanent 
electric dipole moments. The J 2 , V, and W interaction 
terms in Eq. (1) can be understood as the classical dipolc- 
dipole interactions between these permanent dipole mo- 
ments. The J_l interaction term arises due to the transi- 
tion dipole moment between |mo) and \mi). Large chem- 
ical reaction rates [39-41] and large interactions between 
molecules on the same lattice site enforce the hardcore 
constraint [30]. 

The amplitudes and signs of J±, J z , V, and W can 
be tuned independently [30] by tuning the external dc 
electric field and applying external microwave fields [17- 
19, 21-25, 29, 42-46]. The tunneling amplitude t (as- 
sumed to be positive throughout the paper) can be tuned 
by adjusting the depth of the optical lattice. 



II. THE t-J-V-W HAMILTONIAN 



III. ANALYSIS OF THE t-J± HAMILTONIAN 



In the following calculations, we consider nearest- 
neighbor and next-nearest-neighbor interactions on a 
square lattice as shown in Fig. 1 and approximate the 



The simplest experimental realization of Eq. (1) can 
be obtained by applying a very weak external dc electric 
field. In this case, the permanent electric dipole moments 
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TABLE I. (Color online) Symmetries corresponding to the 
irreducible representations of the group D4 as is discussed 
in Appendix A 2 [33, 38, 47]. Wave function symmetries are 
plotted in the xy plane. Note that the A\ representation is 
symmetric under all symmetry operations of the square and 
that the A2 representation is antisymmetric under tv rotations 
about the x and y axes and the lines y = x and y = —x. The 
A2 wave functions can be thought of as being positive (blue) 
on the front and negative (red) on the back, while all other 
wave functions have the same polarity on both sides of the 
wave function. 




0.4 



0.2 



10 



10 



J±/t 



FIG. 2. (Color online) (a): Binding energy A t for two holes 
and other relevant energies Ax, A31, A40, and A2 9 , all of 
which must be positive for the Hamiltonian to support d-wave 

superfluidity within our treatement. (b): A 2^>j for a 
single plaquette, which must be nonzero for the Hamiltonian 
to support d-wave superfluidity within our treatment. 



are very small making J z , V, and W negligible relative to 
J_l, which is proportional to the square of the transition 
dipole moment. In this section, we study the resulting 
t-J± Hamiltonian given by Eq. (1) with J z = V = W = 
and J' z = V = W = 0.' In Sec. Ill A, we describe the 
exact diagonalization of the t-J± Hamiltonian for a sin- 
gle plaquette and identify a set of conditions necessary, 
within our perturbative analysis, for the observation of 
d-wave superfluidity. In Sec. MB, we calculate the phase 
diagram for a two-dimensional lattice of plaquettes per- 
turbatively. 

Throughout the remainder of this paper, we use the 
following notation for states. A ket with one number 
\n) refers to a single plaquette with n = + to- 
tal molecules. A ket with two numbers separated by a 
comma \n^,n^} refers to a single plaquette with spin 
up molecules and spin down molecules. A tensor prod- 
uct of two kets jrifl) \riRi) refers to a system of two pla- 
quettes with tir total molecules on plaquette R and nm 
total molecules on plaquette R' . 



A. Single-Plaquette Analysis 

We use the point symmetries of the square, described 
by the group D4, and the conservation laws of the 
Hamiltonian to simplify the task of diagonalizing the 
single-plaquette Hamiltonian with Hilbert space dimen- 
sion 3 4 = 81 and to understand the symmetries of the re- 
sulting eigenstates. First, the operators and (which 
measure the number of up and down molecules on a sin- 
gle plaquette) commute with the Hamiltonian and with 
each other so we can diagonalize subspaces with fixed 
values of and separately. This reduces the sub- 
space dimensions to 12 for the largest subspaces. We 
then use the basis functions of the group -D4 to diagonal- 
ize the Hamiltonian for each irreducible representation 
separately. This requires diagonalizing 3x3 matrices at 
worst. The symmetries of the resulting eigenstates cor- 
respond to the symmetries of the irreducible representa- 
tions summarized in Tab. I. 

The Hamiltonians for subspaces of constant and 
were diagonalized in the representation basis using the 
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methods described below in Appendix Al. The ground 
states of a single plaquette with fixed total number of 
molecules n = + are summarized in Tab. II. In 
order to construct robust d-wave superfluids at 3/4 filling 
(three molecules per plaquette), we would like it to be 
energetically favorable for two holes to condense on the 
same plaquette. This condition is achieved if the binding 
energy of two holes [33-38] 



A t = 2E g (3)-E g (4)-E g (2) 



(2) 



is positive, where E g (n) is the energy of the ground state 
of \n). A t is shown in Fig. 2(a). The d-wave matrix 
element between |2) and |4) 



must also be nonzero for d-wave superfluidity to be pos- 
sible, is the ci-wave symmetric pair creation operator 
defined as [37, 481 



Jx/t 


Tit) "4. 


Representation 


G (— oo, oo) 


1,0 


Ax 


< -1.22 


1,1 


E 


> -1.22 


1,1 


Ai 


< 


2,1 


A 2 


G (0,0.66) 


3,0 


A 2 


G (0.66,4.62) 


2,1 


E 


> 4.62 


2,1 


Si 


< 


2,2 


Si 


> 


2,2 


Si 



TABLE II. Ground-state configurations and symmetries of 
states with constant n — rif + n±. The results are symmetric 
under interchange of up spins with down spins. The ground 
state for four molecules is always |2, 2) in the Bi representa- 
tion; however, there are two copies of the B\ representation 
for the 1 2, 2) subspace, see Appendix A, and the particular 
ground state changes at J±/t = even though the symmetry 
of the ground state remains d x 2_ y 2. 
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A 



is shown in Fig. 2(b). 



From the information in Tab. II, we can understand 

shown in Fig. 2. We 



the behavior of A t and 



A 



see that for all values of Jx/t, the ground state of |4) 
is a 1 2 , 2) in the B\ representation and thus always has 
d x 2_ y 2 symmetry. For J±/t < —1.22, the ground state 
of 1 2) is a |1,1) in the E representation and thus can- 
not exhibit d-wave superfluidity. For Jx/t > —1.22 the 
ground state changes to a |1, 1) in the A\ representation 
and thus has s-wave symmetry. This change in repre- 



sentations accounts for 



A 



becoming nonzero at 



Jj_/t = —1.22 and the kink in A t at the same value. The 
change in the |4) ground state at J±/t = 0, while retain- 
ing the d x 2_ y 2 symmetry, accounts for the large jump in 



at this value. 

For < Jx/t < 0.66, the ground states of |3) are a 
|3,0) and a |0, 3) in the A2 representation; for all other 
values of Jx/t, the ground states are a |2, 1) and a |1, 2). 
At J±/t = 0.66 the state changes to the E representation 
and at J±/t = 4.62 the state changes to the B\ repre- 
sentation. These changes in representations account for 
the kinks in A t at these values of Jx/t- A f crosses the 
positive J_i/t axis at Jx/t = 0.82. For J ± /t > 0.82, 



the condition ^4 A^ 2j ^ also holds. Thus in this 

regime, holes bind into d-wave symmetric pairs, which 
are necessary for d-wave superfluidity. 

The t-J± Hamiltonian is an improvement over the 
Hubbard Hamiltonian in that the binding energy A t for 
the Hubbard Hamiltonian is positive only for a narrow 
parameter region < U/t < 4.6 [37]. Furthermore, 
the binding energy reaches a maximum for the Hubbard 
Hamiltonian [37] while it increases with J±/t for the t-J± 
Hamiltonian. 

For positive values, the binding energy A t is the 
amount that the ground state |2) |4) or |4) |2) is energet- 
ically favorable over the ground state |3) |3), which could 
be coupled to |2) |4) and |4) |2) through tunneling t' . For 
positive values of A t it is also necessary to consider the 
energy difference between the lowest energy states cou- 
pled to 1 2) 1 4) and |4) |2) through the spin interaction J' ± : 

Ai = E g (0, 2) + E g (3, 1) - E g (l, 1) - E g {2, 2), 

where E g (n^, n^) is the energy of the lowest n^) state. 
Aj^ is shown in Fig. 2(a) for values of J±/t where A t > 

and ^4 A^ 2\ ^ 0. The size of the binding energies 

At and Ai roughly correspond to how large t' and J' ± 
can be, respectively. 

For J±/t > —1.22, the ground state of |2) is an s-wave 
1 1,1) in the Ai representation. However, another |1,1) 
state in the E representation becomes close in energy to 
the ground state for large Jx/t- So we define a third 
energy difference 

A 2g = E g (l, 1(E)) -^(1,1(^0) 

to quantify the energy gap between these two |2) states. 
Here E g (l,l(T)) is the energy of the lowest |1,1) state 
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in the irreducible representation T. A2 g roughly cor- 
responds to how large the overall perturbing Hamilto- 
nian H c g linking plaquettes can be. A 2g is also shown 
in Fig. 2(a) for values of J±/t where A t > and 

At | 2)^0. 

Provided that A t > 0, in a full lattice of plaquettes 
with 3/8 of the lattice sites occupied by up molecules 
and 3/8 of the lattice sites occupied by down molecules, 
it could, in principle, be energetically favorable for the 
molecules to arrange themselves in configurations other 
than two and four molecules on a plaquette throughout 
the entire lattice. Since A t > 0, it will be energetically 
costly to change |2) |4) to |3) |3). Two four molecule pla- 
quettes 1 4) 1 4) cannot rearrange their molecules since they 
are at maximum filling. However, it is possible for two 
two molecule plaquettes |2) |2) to rearrange themselves 
to |3) |1) or |4) |0). Thus it is also necessary to consider 
the binding energies 

A 31 =E g {3)+E g (l)-2E g (2) 



and 



A 40 = E g (4) + E g (0) - 2E g (2) = E g (A) - 2E g {2). 

A31 and A4 are shown in Fig. 2(a). 

As is shown in Fig. 2(a), when A t is positive, the 
other relevant energies are also positive for J±/t less 
than about 12. For J±/t greater than about 12, A40 
becomes negative, and it becomes energetically favorable 
for a 1 2) 1 2) to change to a |4) |0) or |0) |4), in which case 
the manifold of states consisting of only |2) and |4) pla- 
quettes stops being the true ground state. However, the 
corresponding phase diagram can still be studied exper- 
imentally by adiabatically preparing these - no longer 
ground - states beginning with easily preparable excited 
states. A calculation to minimize the energy of a full 
lattice of plaquettes with t' = J' ± = confirms that a 
lattice with half of the plaquettes as |2) and half as |4) is 
the ground state when A40 and A t are positive. 

Although it is outside the scope of the present paper, 
we point out that the narrow region < J±/t < 0.66, 
where |0, 3) and |3,0) are the degenerate ground states, 
could support Nagaoka ferromagnetism [49, 50]. 



B. Double-Plaquette Analysis 

In this section, we describe the behavior of the full lat- 
tice of plaquettes. We use Schrieffer- Wolff transforma- 
tions to find the interactions between nearest-neighbor 
and next-nearest-neighbor plaquettes to second order in 
t' /t and J'j_/t. We then derive an effective XXZ Hamil- 
tonian and solve for the phase diagram in the perturba- 
tive limit. 



1. Schrieffer- Wolff Transformations 

We solve the problem of two coupled plaquettes with 
second-order perturbation theory through the Schrieffer- 
Wolff transformation [51]. There are six relevant doublc- 
plaquette problems in describing a full two-dimensional 
lattice. A |4) can interact with another |4) or with a |2), 
and a |2) can also interact with another |2). For each of 
these three cases, the plaquettes can be situated either 
next to each other or diagonally. 

In all cases, the unperturbed Hamiltonian is the tensor 
product of two single-plaquette Hamiltonians 



H 




U 



(5) 



written here in the basis that diagonalizes it. If the di- 
mension of the two-plaqucttc Hilbert space is d and the 
ground state is Z-fold degenerate, then the upper left is 
an I X I zero matrix and U is a (d — I) x (d — I) diagonal 
matrix with the energy differences between the ground 
and excited states along the diagonal. The |2) and |4) 
energies are different; however, since the number of |2) 
and 1 4) plaquettes is constant, we drop these energies 
here. 

The perturbing Hamiltonian is 



Hi 



H\g Hl ge 



where H\ g defines the first-order shifts in the energies of 
the ground states due to the perturbation, H\ ge defines 
the couplings between the ground and excited states, and 
H\ e defines the couplings between the excited states and 
the first-order shifts in the energies of the excited states. 
The effective Hamiltonian for the low-energy subspace is 
then 



Hcfi — H\ g — HigJJ H lge + 



(6) 



For the remainder of the paper, we divide the Hamilto- 
nians by t. Then the first term is of order t' /t and J' ± /t, 
while the second term is of order (t' /t) 2 and (J' ± /t) 2 . 

We then have the following effective Hamiltonians 
when the two plaquettes are situated next to each other 



rr(4,2) 
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Here the functions / and g on the right hand sides de- 
pend on the interaction strength Jj_/t and describe the 
perturbative coupling of the plaquettes. The superscript 
(n R , n R >) refers to the number of molecules on neighbor- 
ing plaquettes R and R' . The subscripts t and _L refer 
to interplaquette couplings driven by t' and J^_, respec- 
tively. For the t-J±_ Hamiltonian, there are no first-order 
shifts in the ground-state energies so there are no terms 
proportional to J' ± /t. The |2) |4) and |4) |2) states are 
only coupled through tunneling to second order, so there 
is no gf' 2 ^ function. The |4) |4) state cannot couple to 
itself through tunneling since both plaquettes are fully 
occupied so there is no /j 4 ' 4 '' function. 

The effective Hamiltonians for two plaquettes situated 
diagonally are 



H, 



(4,2) 



H. 



cff 
(2,2) 



^(4,2) 


^(2,2) 





^(4,2) 



basis 



12) |4) 
|4) |2> 



off 



basis 1 2) |2) , 
basis 1 4) |4) , 



where 
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h (2,2) = 
^(4,4) = 
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(9a) 

(9b) 
(9c) 



(10a) 
(10b) 
(10c) 



When the plaquettes arc situated diagonally, there is 
no tunneling between them, so there are no terms pro- 
portional to (t'/t) 2 . In particular, the |2) |4) and |4) |2) 
states cannot couple to each other through tunneling in 
second order, so the off-diagonal terms are zero. 



2. XXZ Effective Hamiltonian 

The full lattice of plaquettes can be mapped to an 
XXZ spin Hamiltonian [33, 35, 37, 38] where each pla- 
quette becomes a lattice site, labelled by R, and the 
states 1 2) and |4) of two and four molecules become the ef- 
fective spin up and spin down states, respectively. 



Using the functions Eq. (8) and Eq. (10) calculated using 
the Schrieffer- Wolff transformation, the effective Hamil- 
tonian is 

{R,R') 

5 (4 ' 2) (S+Sr, + S R S+)+f^n Ri[ n m +f^n m n R% 



i(R,R')) 



+ hS 2 ' 2) n m n R ^ + h^nm.riR'u. 



Since each site R has either one spin up or one spin down, 



1 



n Rf = 2 + S R and 
Thus, dropping constant terms, 



flRJ). 



1 -S> 



H cfi — ^2 \j± (SrSri + SrSr') + J z iS R S 



{R,R') 



J Z 2S z R S z R , + BY, S Ri (ii) 

((R,R'» R 



where 



J±=2g^ 2 \ (12a) 
J,i = / (2 ' 2) +/ (4 ' 4) -2/^, (12b) 
J z2 = hW +h&*> ~2h^ 2 \ (12c) 

B = 2 (,/( 2 ' 2 ) - + - . (12d) 

Since we are interested in the phase diagram at constant 
3/4 filling of the lattice with molecules, J2r^r — is 
constant in the XXZ effective Hamiltonian, so we can 
neglect the magnetic field term B. In fact, for the same 
reason, we have already dropped the energies of |2) and 
|4) in H [Eq. (5)]. 

3. Phase Diagram for the t-J± Hamiltonian 

If we only consider nearest-neighbor interactions, then 
Jz2 = 0. In this case, there are three phases, which 
can be qualitatively understood by considering the fol- 
lowing limits of the XXZ model [52]. If J zl > |Jj_|, 
it is energetically favorable for the spins to be anti- 
ferromagnetically ordered in the Z direction correspond- 
ing to a checkerboard solid of bound holes. If J z i <C 
— |jj_|, it is energetically favorable for the spins to be 
ferromagnetically ordered in the Z direction; however, 
since J2 R $r i s fixed, this corresponds to a phase sepa- 
ration of the spins and a phase separation of the bound 



holes. Finally, if I Jj_| > 



Jz 



it is energetically favor- 



able for the spins to be ordered in the XY plane, which 
corresponds to a superfluid of bound holes. Specifically, 
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FIG. 3. (Color online) Phase diagram for the t-J± Hamilto- 
nian. The region to the right of the vertical black line is the 
region where A t > 0. The line J'±/t' = J±/t is shown as a 
dashed green line and passes through all three phases. 



for Jj_ < (Ji > 0), the spin order is ferromagnetic 
(anti- ferromagnetic) , corresponding to a d-wave super- 
fluid with correlation function (AL^A^h') whose sign 
has uniform (checkerboard) structure in the R—R' plane. 
Since Jj_ can be mapped to — J±_ by a sublattice rotation 



Jzl 



between the 



[52], we do not distinguish, for | Jj_| ^> 

ferromagnetic and anti-ferromagnetic cases and simply 
refer to both phases as a d-wave superfluid. The phase 



transitions occur at 



J 



Jzl 



[52, 53]. 



The case of nonzero next-nearest-neighbor interac- 
tions, nonzero J z2 , has been studied numerically in 
Ref. [53] and with mean field theory in Refs. [54-56]. 
For example, in Ref. [53], it is shown that, for a cer- 
tain parameter range satisfying J z2 > J z i > and 
Jz2 J± > 0, it is energetically favorable for the pla- 
quettes to arrange themselves in a striped solid. Assum- 
ing, by analogy, that \J z2 \ > \J z i\ and \J z2 \ > \J±\ are 
both necessary for a new phase to appear, no such phase 
can occur in the perturbative phase diagram for the t- J± 
Hamiltonian since this set of conditions is never satis- 
fied. Furthermore, near the phase transition boundaries 
for this phase diagram, the magnitude of J z2 is about an 
order of magnitude smaller than the magnitude of J z \. 
We therefore expect J z2 to have an insignificant effect on 
the locations of phase transitions, so we neglect J z2 for 
the remainder of the paper. 



The condition for a phase transition ± 
thus 



= Jzi is 



Here, a transition between superfluid and checkerboard 
solid occurs for +, a transition between superfluid and 
phase separation occurs for — , and the functions are eval- 
uated at J±/t. By solving Eq. (13), the phase diagram in 
the original variables Jj_ jt and J' ± jt' is computed in the 
region At > and is shown in Fig. 3. All three phases 
are present in this phase diagram and the d-wave super- 
fluid phase exists for a large range of values of J± jt and 
J'j_/t' . The easiest case to study experimentally is t! = t 
and J'x = J±, which is outside the validity of this pertur- 
bative calculation. However, as a guess as to the physics 
for these values of t' and J' ± , it is useful to consider the 
line J'j_/t' = J±/t. As is shown in Fig. 3, this line passes 
through all three phases, indicating that all three phases 
might be observable in the simplest t- J± experiment with 
a homogeneous lattice. 

There is a strong indication that the qualitative 
features of our results may be relevant to the non- 
perturbative regime where J'^/tf = J±/t. Specifically, 
the phase diagram along the line J'j_/t' = Jj_/t in Fig. 3 
is qualitatively similar to the phase diagram along the 
line of 1/4 hole density in Fig. 4 of Ref. [57], which nu- 
merically studies the standard t-J model. Indeed, the 
order of the phases in the two diagrams is the same pro- 
vided that one identifies the region of uncondensed bound 
holes in Ref. [57] with our checkerboard solid phase and 
the Fermi liquid in Ref. [57] with our region A f < 0. 



IV. ANALYSIS OF THE EFFECTS OF J z , V, 
AND W 

In this section, we analyze the effects of J z , V, and 
W on the d-wave superfluid phase. First in Sec. IV A, 
we repeat the single-plaquette analysis of Sec. Ill A for 
various Hamiltonians with nonzero J z , V, and W. Then 
in Sec. IV C, we compute the phase diagrams using the 
methods of Sec. IIIB for those Hamiltonians capable of 
exhibiting d-wave superfluidity 



A. Single-Plaquette Solutions 

First, we examine the Hamiltonian with independent 
Jj_ and J z and zero V and W, the t-J±-J z Hamiltonian. 
Contour plots for the binding energy A t and {^4 A^ 2^ 
for this Hamiltonian are shown in Fig. 4. From these 
figures we see that A t > and ^4 A^ 2^ 7^ for most 

values of J±/t > and J z /t > 0. 

The 1 4) ground states in the A 2 and B\ representations 
for positive J z /t are 



±2 



9i 



(4,2) 



fi 



(2,2) 
t 

v 



(4,2) 



(2,2) 



- 2/r ; 



(13) 



|2,2(A a )) = -^(itm>-imt>), 

12, 2(Bx)> = (itm) + lint)) + o (?f\ 
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FIG. 4. (Color online) Single-plaquette analysis for the t- 
J±-J z Hamiltonian. (a): Binding energy A t . Contours of 
constant At/t are labeled, (b): Matrix element |^4|a^| 2^|. 

Contours of constant (4 At 2) are labeled. 



FIG. 5. (Color online) Single-plaquette analysis for the the t- 
J-W Hamiltonian with J± = J z = J. (a): Binding energy At. 
Contours of constant At/t are labeled, (b): Matrix element 

(4 2) . Contours of constant (4 Aj 2) are labeled. 



Thus, for positive J z /t and along the line J±/t = 0, the 
two alternating spin configurations are the two-fold de- 
generate 1 4) ground states. When J±/t is made nonzero 
but small, the symmetric combination of these two states, 
which has B\ symmetry, becomes the non-degenerate 
ground state. 

Next, we examine the Hamiltonian with independent 
J± = J z — J and W and zero V, the t-J-W Hamiltonian. 
Contour plots for the binding energy A t and ^4 A^ 2^ 
for this Hamiltonian are shown in Fig. 5. From these 



figures we see that A t > and ^4 A^ 2j =^ for most 

values of J/t > and W/t approximately between the 
lines W/t = ±2 J/3t. 

This behavior can be explained by considering the 
ground-state configurations and symmetries of \2) and 
1 4) shown in Fig. 6. If W is large and positive, all of 
the spins will point down in the single-plaquette ground 
states: the \2) ground state is a p- wave |0, 2) and the |4) 
ground state is a ci-wave |0, 4). For W = and J/t > 0, it 
is energetically favorable to have an equal number of spin 
up and spin down molecules on each plaquette, and the 
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FIG. 6. Configurations of the single-plaquette ground states 
for two and four molecules for the t-J-W Hamiltonian with 
J± = J, = J. 



|2) ground state is an s-wave |1, 1), while the |4) ground 
state is a d-wave |2,2). As W is decreased from a large 
positive value (and J±/t > 0), the |2) ground state first 
switches from |0, 2) to an s-wave |1,1), and then the |4) 
ground state switches from |0,4) to an s-wave |1, 3). As 
W is further decreased, the |4) ground state switches to 
a d-wave symmetric \2, 2). As W is made large and neg- 
ative, the process repeats with up spins replacing down 

spins. ^4 2^ vanishes for W/t approximately out- 
side the lines W/t = ±2J/3t since the |4) ground state 
switches from d-wave to s-wave outside this region. In 

particular, ^4 2^ = in the narrow region where the 
ground states are s-wave |1, 1) and d-wave |1,3) because 
A t creates one up spin and one down spin. 



We do not consider Hamiltonians with nonzero V here 
since they do not support d-wave superfluidity in our 
perturbative calculations as is described in the following 
section. We also note that, similar to the t-J± model 
where A40 becomes negative at large J±/t, the region 
where a lattice of |2)'s and |4)'s is the true ground state 
is also limited in the models discussed in this section. 
However, as in the case of the t- J± model, the full phase 
diagrams discussed in Sec. IV C can still be accessed in 
experiments with ultracold polar molecules by adiabatic 
preparation starting with suitable initial states. 



B. First-Order Contributions to the XXZ Effective 
Hamiltonian 

In this section, we examine the question of which of the 
Hamiltonians defined by Eq. (1) have nonzero first-order 
contributions to the XXZ effective Hamiltonian. A first- 
order shift in the ground-state energies, a nonzero H\ g 
in Eq. (6), would contribute to, and generally dominate, 
the diagonal terms f^ 2 \ f^ 2 \ f^ 4 \ h^ 2 \ h^ 2 \ and 
/j(4,4) no t the off-diagonal term g^ 4 ' 2 \ Since the or- 
der (f /t) 2 function </ 4,2 ' is the only function contribut- 
ing to Jj_ [Eq. (12a)], J± will generally be small com- 
pared to J z \ and J Z 2 [Eqs. (12b) and (12c)], if there 
are nonzero first-order contributions. Thus, this per- 
turbative analysis predicts no robust superfluid phase if 
there are nonzero first-order contributions to the effec- 
tive Hamiltonian. Therefore, we will not compute the 
phase diagram for those Hamiltonians that contribute 
first-order corrections. Note that there are no first-order 
shifts in the t-J± Hamiltonian studied in Sec. III. 

To study when first-order contributions arise, we con- 
sider here the effects of the J' z , V, and W' terms sepa- 
rately. Since these terms cannot couple |2) |4) to |4) |2), 
they will only contribute to the diagonal matrix elements 
of fleff- Let R and R' label the plaquettes containing 
sites r and r' , respectively, and let n Ra denote the num- 
ber of molecules with spin a on plaquette R. Due to the 
D4 symmetry, any single-plaquette eigenstate I^Rf, 
that is non-degenerate within the manifold of states with 
constant nn\ and n R ^ satisfies 

(n R ^,n Ri \n ra \ n m , n Rl ) = ^n Ra . (14) 

First, consider the V' term. From Eq. (14), the state 
I^-r) \tir>) satisfies 

(n r n r ,) = ~~Ti R n R / , 

where n r = + n r ^ and n R = n R -f + n R ^. Thus the V' 
contributions vanish to first-order if and only if at least 
one of the two interacting plaquettes is empty. Since the 
effective Hamiltonian is constructed from |2) and |4) pla- 
quettes, the V term will always contribute to first order 
and will not be considered further here. It is important 
to emphasize that our analysis should be considered to 
be exactly valid only in our perturbative regime, since it 
is believed that the t-J model, which contains nonzero 
V, supports d-wave superfluidity [1-3, 57, 58]. 

Next, consider the J' z term. Since S* — (n r f — n r ±)/2, 
the state \n R ^,n Ri _) \n R ^,n R ^) satisfies 



(s;sz 



1 

04 



(n Rt - riRi) (riH't - n R'l) 



Thus, the J' z contributions vanish to first-order if and 
only if at least one of the two interacting plaquettes sat- 
isfies n-f = nj,. In all of the regions identified in Sec. IV A 
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as possibly supporting d-wave superfluidity in the t-J±- 
J z Hamiltonian, \2) and |4) have this property. There- 
fore, J' z never contributes at first-order in the parameter 
regimes that we are interested in for this Hamiltonian. 

Finally, consider the W term. The state 
\n R ^,n Ri ) \n R , t ,n Rli ) satisfies 

(n r S*, + n r ,S z r ) = -^{n m n R ^ - nR±n R i\). 

Thus the W contributions vanish to first-order if and 
only if the two interacting plaquettes satisfy n R ^n R ^ — 
n R l n R 'l- From Fig. 6, we see that this condition is met 
for 1 2) 1 4), 1 2) 1 2) and |4) |4) only for the region where the 
|2) ground state is s-wave symmetric |1,1) and the |4) 
ground state is d-wave symmetric |2,2). We also note 
that the condition = n^, necessary for the J' z contri- 
butions to vanish at first-order, is satisfied by both |2) 
and |4) in this region. Notice, however, that the condi- 
tion = ni is not satisfied by the |4) states outside this 
region. Thus, not only W but also J' z will give nonzero 
first-order contributions outside this region. 

While other regions in Fig. 6 cannot exhibit a d-wave 
superfluid within our analysis, they may still exhibit in- 
teresting phases at appropriate filling fractions. For ex- 
ample, the regions where the ground states are |1, 1) and 
|0,4) or |4,0) might support a d-wave solid phase with 
an asymmetry between up and down spins, while the re- 
gions where the ground states are |1, 1) and |1, 3) or |3, 1) 
might support an s-wave solid phase. However, we will 
not discuss such phases further and will focus, instead, 
on the parameter space that has no first-order correc- 
tions and that is therefore capable of exhibiting d-wave 
superfluidity. 



C. Phase Diagrams 

We first consider the general t-J±-J z Hamiltonian. 
There are four independent parameters, J±/t, J z /t, 
Jj_/t ', and J' z /t', so the phase diagram is four- 
dimensional. Notice that, within the perturbative treat- 
ment, the fifth parameter t' /t affects only the overall en- 
ergy scale but not the phase diagram. In order to plot 
a manageable phase diagram, we restrict the parameter 
space to J'j_/t' = J±/t and J' z /t' = J z /t. Using the ex- 
pressions for Jj_ and J z \ given by Eq. (12) and the gener- 
alization of Eq. (8) to include J z and J' z , we find the phase 



boundaries by solving 



J± 



The resulting phase 



diagram is shown in Fig. 7. As is discussed in Sec. IV A, 
along the J z /t axis, the |4) ground state becomes dou- 
bly degenerate. Thus the effective Hamiltonian no longer 
maps to an XX Z Hamiltonian, and we do not compute 
the phase diagram along this line. In the absence of the 
hardcore constraint, the t-J z Hamiltonian on a homoge- 
neous lattice is studied in Ref. [32]. The t-J± phase dia- 
gram along the green dashed line J' ± /t' = J±/t shown in 
Fig. 3 corresponds to the phase diagram along the J±/t 




phase separation 
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checkerboard 
solid 
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FIG. 7. (Color online) Phase diagram for the t-J±-J z Hamil- 
tonian assuming J'^/t' = J±/t and J' z /t' — J z /t. The region 
above the black curve is the region where A t > 0. The SU(2)- 
symmetric Hamiltonian with J± — J z is shown by the dashed 
green line. We do not compute the phase diagram along the 
J z /t axis since the |4) ground state is degenerate here. 



axis in Fig. 7. For values of J z /t greater than about 1, 
the checkerboard solid phase is no longer present, and 
for values of J z /t greater than about 2, the superfluid 
phase is no longer present. We see that large values of 
J z have the effect of reducing the areas of the checker- 
board solid and superfluid phases. However, by tuning 
J z /t to a value between about 1 and 2, the extent of the 
d-wave superfluid phase along the J±/t direction is in- 
creased relative to its value without the J z interaction. 
This is achieved by the suppression of the checkerboard 
solid and phase separation phases at small J±/t and large 
J±/t, respectively. 

We next consider the phase diagram for the SU(2)- 
symmctric t-J±-J z Hamiltonian with J± — J z = J shown 
in Fig. 8 as solid lines. The phase diagram along the green 
dashed line J'/t' = J /t in Fig. 8 for the SU(2)-symmetric 
t-J±-J z Hamiltonian corresponds to the phase diagram 
along the green dashed line J± — J z in Fig. 7. The 
regions not in phase separation are reduced from those 
for the t-J± Hamiltonian (see Fig. 3). The fact that 
J z reduces the superfluid phase for the SU(2)-symmetric 
case is to be expected from Fig. 7 since the line J± = J z 
does not pass through the regions where the superfluid 
phase is enhanced by the presence of J z . 

Finally, we add W and consider the phase diagram 
for the t-J-W Hamiltonian with J± = J z = J shown 
in Fig. 8. The phase diagram is shown along the lines 
W = (discussed above), W = J/4, and W = J/2, 
which are contained within the region where |2) is an s- 
wave symmetric |1, 1) and |4) is a d-wave symmetric |2, 2) 
(see Fig. 6). Increasing W moves the transition between 
the superfluid and phase separation phases up along the 
J'/t' axis and therefore slightly increases the region of 
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FIG. 8. (Color online) Phase diagram for the t-J-W Hamil- 
tonian for J± = J z = J . Phase transitions for W = (the 
SU(2)-symmetric t-J±-J z Hamiltonian) are shown in solid 
lines, transitions for W = J/4 are shown in dashed lines, 
and transitions for W = J/2 are shown in dotted lines. The 
regions to the right of the vertical black lines are the re- 
gions where At > for these three values of W. The line 
J' /t' — J/t is shown as a dashed green line. 

d-wavc superfluidity. Increasing W also decreases the 
extent of the checkerboard solid phase along the J/t axis. 

In summary, within our treatment, Hamiltonians in- 
volving V do not support superfluidity. At the same 
time, we identify regions of parameter space where J z 
and W enhance the d-wave superfluid phase. 

V. EXPERIMENTAL REALIZATIONS OF THE 
PERTURB ATIVE CALCULATION 

The perturbative results obtained above can be re- 
garded only as a qualitative guess as to the behavior of 
the simplest homogeneous square lattice since intrapla- 
quette and interplaquette couplings are equal in this 
case. Therefore, in this section, we propose experimen- 
tal configurations accurately described by our perturba- 
tive analysis. In Sec. V A, we describe a one-dimensional 
stack of plaquettes. Then, in Sec. VB, we briefly describe 
an experimentally more challenging two-dimensional con- 
figuration. 

A. One-Dimensional Stack of Plaquettes 

As shown in Fig. 9, we propose to apply a dc electric 
field in the vertical (z) direction. We then propose to 
stack plaquettes with sides of length a on top of each 
other along z a distance a' apart. The strength of dipolar 
interactions between two molecules is proportional to 




J,t 

FIG. 9. (Color online) Geometry of the one-dimensional stack 
of plaquettes. The external dc electric field is perpendicular 
to the plaquettes. 



where r is the distance between the molecules and 9 is 
the angle made between the line connecting the molecules 
and the external dc electric field. Due to the 1/r 3 depen- 
dence, interplaquette dipolar interactions can be treated 
perturbatively relative to the intraplaquette interactions 
provided a' 3> a. One way of achieving this geometry 
experimentally is to use different frequency lasers to cre- 
ate the optical lattices in the plane and in the vertical 
direction. The distances a and a' can also be controlled 
by varying the angle at which the lasers interfere [59, 60] 
or by holographic techniques [61]. 

The plaquettes can be made by interfering lasers of 
wavelength 2a and 4a to create a superlattice [9] such 
that tunneling between plaquettes in a plane is negligi- 
ble. To avoid in-plane interplaquette dipole-dipole inter- 
actions, molecules in neighboring stacks may have to be 
removed. The required addressability can be acheived 
by applying temporary additional light shifts or electric 
field gradients. Alternatively, instead of emptying neigh- 
boring stacks, an extreme version of the superlattice can 
be used to separate the stacks enough to make both the 
tunneling and the dipolar interactions between them neg- 
ligible. 

As is shown in Fig. 9, with this geometry, there is a 
hopping amplitude t between nearest neighbors within 
a plaquette and a hopping amplitude t' along z between 
nearest-neighbor plaquettes. The amplitudes t and t' can 
be controlled separately by varying the intensities of the 
lasers making the lattices in each direction. There are 
five strengths, two within a plaquette and three between 
plaquettes, of the dipolar interactions J±, J z , V, and W 
discussed here generically as J. As shown in Fig. 9, let 
J be the strength of dipolar interactions between nearest 
neighbors within a plaquette, and let J', J", and J'" be 
the strengths of the dipolar interactions between plaque- 
ttes. The ratios J'/ J, J" /J, and J'" /J are functions of 
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FIG. 10. (Color online) Phase diagram for the t-J± Hamil- 
tonian for a one-dimensional stack of plaquettes with a' /a = 
5/2. Regions to the right of the vertical black line are where 
A t > 0. The lines -J'Jt! = J±/2t, -J'Jt' = J±/t, and 
— J'±/t' — 2J±/t are shown by the green dashed lines. 



the ratio a' /a and are controlled separately from the ra- 
tio t'/t. For a' > a, by Eq. (15), J', J", and J'" are the 
opposite sign of J. 



1. Phase Diagrams 

For a one-dimensional stack of plaquettes, we calculate 
the same phase diagrams computed above for the two- 
dimensional lattice of plaquettes. In an experiment with 
one independent dipolar interaction strength J, there are 
three independent parameters: J/t, a' /a, and t'/t or, 
equivalently, J/t, J' /t' ', and a' /a. In order to compare 
with the phase diagrams computed above, we choose the 
latter set of parameters. It is currently possible to use 
lasers of wavelength 1064 nm to produce optical lattices 
with spacing a' = 532 nm. Assuming that it is also possi- 
ble to use a second wavelength in the range 400 — 600 nm, 
in the following phase diagrams, we use a' /a = 5/2. 
For a' /a = 5/2, J'/J = -0.128, J"/ J = -0.081, and 
J'" /J = -0.054. 

The phase diagram for the t-J± Hamiltonian is shown 
in Fig. 10. Since it is possible with this scheme to con- 
trol J'Jt! by varying a! /a and t'/t, we plot the lines 
-J'Jt' = J x /2t, -J'Jt' = J x /t, and -J'Jt' = 2J x /t 
in Fig. 10. This phase diagram is qualitatively similar to 
the corresponding phase diagram for the two-dimensional 
case (Fig. 3). However, the overall vertical scale of the 
phase diagram is increased for the stack resulting in a 
larger region of d-wave superfluidity The increase in the 
vertical scale can be explained by the presence of the J" 
and J'" terms, which reduce the f± functions relative to 
the ft and gt functions in Eq. (8). Thus, a larger value 
of J' ± /t' is needed to reach the phase boundaries. We 
refer the reader to Appendix A 3 for a symmetry-based 



FIG. 11. (Color online) Phase diagram for the t-J±-J z Hamil- 
tonian assuming J'±/t' = J±/t and J' z /t' = J z /t for a one- 
dimensional stack of plaquettes with a /a — 5/2. The region 
above the black curve is the region where At > 0. The SU(2)- 
symmetric Hamiltonian with J± — J z is shown by the dashed 
green line. 



explanation for why J" and J"' reduce the f± functions. 
Similar arguments hold for the increased scales in the 
t-J±-J z and the t-J-W phase diagrams discussed below. 

The phase diagram for the t-J±-J z Hamiltonian with 
J'Jt! = Jj_/t and J'Jt' = Jjt is shown in Fig. 11. This 
phase diagram is qualitatively similar to the correspond- 
ing phase diagram for the two-dimensional case shown in 
Fig. 7. As with the two-dimensional case, J z can increase 
the superfluid phase. 

The phase diagram for the SU(2)-symmetric t-J±-J z 
Hamiltonian is shown in Fig. 12 as solid lines. Again, 
the diagram is qualitatively similar to its counterpart for 
the two-dimensional system (Fig. 8) with the range of 
J' ft' not in the phase separation regime increased. 

The overall scale of the phase diagram in Fig. 11 is 
increased from the corresponding diagram for the two- 
dimensional system shown in Fig. 7 as is expected from 
Figs. 10 and 12. The t-J± phase diagram along the mid- 
dle green dashed line —J' ± /t' = J±/t shown in Fig. 10 
corresponds to the phase diagram along the J±/t axis in 
Fig. 11. The SU(2)-symmetric phase diagram along the 
middle green dashed line —J'/t' — J/t shown in Fig. 12 
corresponds to the phase diagram along the green dashed 
line J_l = J z in Fig. 11. Since both of these green lines 
in Figs. 10 and 12 pass through larger regions of the su- 
perfluid phase than in the corresponding phase diagrams 
for the two dimensional plane (Figs. 3 and 8), the over- 
all scale of the t-J±-J z phase diagram (Fig. 11) increases. 
The superfluid phase can be further increased by decreas- 
ing the slope of the line relating J'j_/t' {J'/t') and J±/t 
(J/t) in Fig. 10 (Fig. 12) since lines with shallower slopes 
pass through larger regions of the superfluid phase. 

The phase diagram for the t-J-W Hamiltonian for 
J± = J z = J is shown in Fig. 12. In contrast to the 
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FIG. 12. (Color online) Phase diagram for the t-J-W Hamil- 
tonian for J± — J z — J for a one-dimensional stack of pla- 
quettes with a' /a = 5/2. Phase transitions for W=0 (the 
SU(2)-symmetric t-J±-J z Hamiltonian) are shown in solid 
lines, transitions for W = J/4 are shown in dashed lines, 
and transitions for W = J/2 are shown in dotted lines. The 
regions to the right of the vertical black lines are the re- 
gions where At > for these three values of W. The lines 
-3' It' = J/2t, -J' /t' = J/t, and -J'/t' = 2J/i are shown 
by the green dashed lines. 



phase diagram for the two-dimensional t-J-W Hamilto- 
nian, in the stack geometry, an increase in W moves the 
transition between the superfluid and phase separation 
phases down rather than up, thus reducing the superfluid 
phase. Furthermore, in the stack geometry, an increase 
in W only slightly decreases the extent of the checker- 
board solid phase. Thus for the stack of plaquettes, the 
W term can slightly increase the superfluid region only 
for small \J'/t'\. 



2. Preparation and Detection 

Since there are no strong relaxation mechanisms in op- 
tical lattice experiments using cold atoms and molecules, 
it is not an easy task to prepare the ground state of a 
given Hamiltonian. One strategy is to prepare an easier 
state that is the ground state of another Hamiltonian and 
to adiabatically change that Hamiltonian to the desired 
one [9, 48, 62, 63]. Here we consider the adiabatic prepa- 
ration of the d-wave superfluid phase at the point where 
J z i = for the t-J± Hamiltonian. This corresponds to 
preparing the ground state of the XX magnet. We pro- 
pose first adiabatically preparing |2) and |4) states alter- 
nating in a stack and then adiabatically preparing the 
ground state of the t-J± Hamiltonian. The preparation 
of the ground state of the XX magnet is similar to the 
method described in [62] for adiabatically preparing the 
ground state of the Heisenberg antiferromagnet. 

First, consider bringing the plaquettes far apart so as 



to avoid interplaquette interactions and preparing the |2) 
and 1 4) ground states alternating along the stack. Let 
H be the t-J± Hamiltonian on a single plaquette and 
consider the following single-plaquette Hamiltonians 

H 2 (t) = H -B(t)(S?-SI), (16) 
H 4 (t) = H - B(t) (SI + S| - SI - St) , (17) 

used to prepare the |2) and |4) ground states, respec- 
tively Here B(t) is the strength of an effective alter- 
nating magnetic field on the sites of the plaquette as a 
function of time r. For large positive values of B(r), the 
ground state of H 2 for two molecules is up on 1 and down 
on 3, while the ground state of H4 for four molecules 
is up on 1 and 2 and down on 3 and 4. These states 
can be prepared using single-site addressability provided 
by electric field gradients or high-resolution optics [6, 7]. 
For values of J±/t of interest, in a plaquette with two 
molecules, the energy gap between the ground state and 
first excited state of H 2 never closes as B(t) is reduced 
to zero. The same is true of H± for a plaquette with four 
molecules. Thus the |2) and |4) ground states of the t-J± 
Hamiltonian can be prepared by adiabatically reducing 
the effective magnetic field B(t) from some large initial 
value to zero in H 2 and f/4, respectively. Once the |2) 
and |4) ground states have been prepared, the plaque- 
ttes can be brought to a distance a' from each other by 
changing the angle at which the lasers interfere or by 
holographic techniques. 

Once the |2) and |4) ground states have been prepared, 
the ground state of the t-J± Hamiltonian at the point 
Jzi — can be prepared by adiabatically turning off the 
effective magnetic field B(t) acting on the effective spins 
I "ft) and IJL) in the Hamiltonian 



H(t) 



J_l {SrSr+i + SrS 



+ B(t)J2(-1) R Sr- (18) 

R 

The first term is just the XXZ effective Hamiltonian 
Eq. (11) at the point J z \ = 0, while the second term 
describes an alternating effective magnetic field between 
plaquettes along the Z direction. The energy gap be- 
tween the ground state and the first excited state of the 
Hamiltonian Eq. (18) decreases monotonically with de- 
creasing magnetic field but never closes, so the d-wave 
superfluid phase can, in principle, be prepared by adi- 
abatically turning off B(t). However, for a stack of TV 
plaquettes, the energy gap between the ground state and 
the first excited state of the XX chain is proportional to 
1/7V. Although it is outside the scope of the present pa- 
per, the ideal B(t) over a given time T can be calculated 
by mapping Eq. (18) to free fermions and maximizing the 
superfluid order parameter when B(T) = 0. Since the en- 
ergy gap is large for large B[t) and only scales as 1/N 
for small B(t), it is likely that this ideal B(t) decreases 
rapidly initially and more slowly at later times. 
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As a rough estimate for the minimum time T necessary 
to prepare the superfluid phase at the point J z \ = 0, we 
estimate T — N/Jx- The strength of the dipole-dipole 
interactions at a distance of 200 nm is roughly Jj_ ~ 
2tt x 4 kHz for the molecule KRb. The value of J ± ~ 
—2ir x 500 Hz is then fixed by the ratio a' /a. t is chosen 
such that J± is as large as possible while perturbation 
theory is valid, t' is chosen such that J z \ — 0. Then 
for J±/t = 4, t ps 2tt x 1 kHz and J z \ = requires 
t' ss 2tt x 300 Hz. At this point, Jx ps 2tt x 1.5 Hz so T ps 
100A^ ms. For these parameters, A t ps 2tt x 2 kHz and 
Aj_ ps 27r x 4 kHz so the conditions t' <C A t and <C 
Aj_ are met, and perturbation theory is valid. Using the 
inverse gap at B(t) = in Eqs. (16) and (17) gives 400 ^s 
as a rough estimate for the minimum time necessary to 
prepare the \2) and |4) ground states. Thus, as expected, 
the preparation time is dominated by the preparation of 
the ground state of the XX chain. Therefore, assuming 
an optimistic coherence time of 1 s, roughly 10 plaquettes 
can be prepared. For LiCs, J± ~ 2tt x 400 kHz at a 
distance of 200 nm allowing roughly 10 3 plaquettes to 
be prepared in 1 s. These numbers can be improved by 
reducing the ratio a' /a or by preparing a superfluid away 
from the point J z \ = 0. 

The effective magnetic fields in Eqs. (16)— (18) can be 
created using tensor shifts [18, 21, 24, 30, 43, 64] and su- 
perlattices [9] so that up and down spins (both the initial 
|t) and \l) and the effective |f|-) and |-|J.)) have different 
potential energies. 

The d-wave superfluid phase can be detected via 
second-order noise correlations in the expanding molec- 
ular cloud which is proportional to the the four-point 
function 

G a cr' (Qi Q') « MQ)rv(Q')) - (n a (Q)) (rv(Q')) 

at the time the molecules are released from the trap [37, 
65, 66]. Here 

n o (Q) oc^e'^ct'V. 

is the quasi-momentum distribution and L rr / is the vec- 
tor connecting lattice sites r and r'. Since G-t^(Q,Q') 

contains terms proportional to ^A^ ^A,^/ ^ (where R 

and R' label plaquettes) , a ci-wave superfluid will exhibit 
interference fringes at Q + Q' = 2itmi / a' for any integer 
(half-integer) m where z is the unit vector along the stack 
when Jx < (> 0). Since the state has d x 2_ y 2 symmetry, 
these fringes will be modulated in the x-y plane by an en- 
velope that vanishes along the nodal lines Q x = ±Q y and 
Q' x = ±Q' y [37]. Therefore, the behavior of the fringes 
will be similar to that shown in Ref. [37], except that 
the superfluid phase will be signaled by fringes along i 
instead of fringes in the x-y plane. Due to the absence 
of cycling transitions in molecules, the noise correlation 
measurements may have to be done by first converting 
the molecules back into atoms [12]. 



B. Two-Dimensional Realization of the 
Perturbative Regime 

While we have discussed the one-dimensional stack 
of plaquettes, d-wave superfluidity is associated with a 
plane. An experiment observing <i-wave superfluidity in 
a two-dimensional lattice would therefore be more rele- 
vant. The simple solution of increasing the spacing be- 
tween plaquettes in a two-dimensional lattice to reduce 
the dipolar interactions is not feasible since this would 
typically reduce the tunneling amplitude t' between pla- 
quettes to essentially zero. 

However, if the energy difference between the rotor lev- 
els is not constant throughout a lattice of plaquettes, then 
the Jx dipolar interaction between plaquettes could be 
suppressed [67]. This can be achieved by introducing an 
extra splitting A between rotor levels in some plaque- 
ttes and not in others. Then the Hamiltonian for each 
plaquette with the extra splitting gets an extra term 

r 

where the sum is taken over the four sites of the pla- 
quette. Let (R x ,R y ) be the integer coordinates of the 
plaquette labeled by R. Suppose that the splitting is ar- 
ranged in a checkerboard fashion, so that for integers n 
and m plaquettes (R x + 2n, R y + 2m) get the extra split- 
ting A while plaquettes (R x + 2n + 1, R y + 2m + 1) do 
not. In this case, dipolar interactions between nearest- 
neighbor plaquettes will be suppressed by A since they 
are off resonant; however, dipolar interactions between 
next-nearest-neighbor (diagonal) plaquettes will not be 
suppressed since they will stay resonant. Thus to make 
the perturbative calculations valid, four separate split- 
tings would need to be introduced to suppress both near- 
est and next-nearest-neighbor interactions. However, 
with only one splitting A, a one-dimensional chain of 
plaquettes can be simulated experimentally. 

Finally, for A that is large enough to make perturba- 
tion theory valid, a lattice of plaquettes filled with |1, 1) 
and 1 2, 2) states will typically be an excited-state of the 
system of decoupled plaquettes, and the d-wave super- 
fluid state will typically be an excited state of the full 
Hamiltonian of weakly coupled plaquettes. Nevertheless, 
such states can still be prepared adiabatically from ap- 
propriate excited states. 



VI. CONCLUSIONS 

We have shown that the t- Jx Hamiltonian on a square 
lattice, in the regime of weakly coupled plaquettes, ex- 
hibits a d-wave superfluid phase in addition to the 
checkerboard solid phase and phase separation. The ad- 
dition of large J z or W interactions destroys the super- 
fluid phase; however, we have identified ranges of these 
parameters, for which the superfluid phase is enhanced. 



Any nonzero V destroys the superfiuid phase in this per- 
turbative analysis. 

These perturbative calculations can be used as a quali- 
tative guess for the behavior of the simplest experiments, 
which are outside of the perturbative limit. Further- 
more, the perturbative regime can be accessed in exper- 
iments in a one-dimensional stack of plaqucttcs. The 
phase diagrams for the one-dimensional stack of pla- 
quettes are qualitatively similar to the phase diagrams 
for the two-dimensional lattice of plaquettes. By ex- 
perimentally observing the phase diagrams in both the 
perturbative and non-perturbative regimes for the one- 
dimensional stack, one may be able to understand the 
relationship between the calculations presented here and 
the non-perturbative phase diagrams. This knowledge 
may then be useful in understanding the relationship 
between the two-dimensional phase diagrams presented 
here and experiments on a homogeneous two-dimensional 
lattice. Similar results might be achievable experimen- 
tally with the large-magnetic-moment atoms dysprosium 
[68] and chromium [69] instead of polar molecules. 

Ultracold polar molecules have the potential for ex- 
perimentally observing <i-wave superfluidity in a con- 
trolled environment that could allow us to learn about 
the physics of the t-J model. We hope that the insight 
gained by these investigations could help to explain the 
physics of high-temperature superconductivity and result 
in many theoretical and practical applications. 
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Appendix A: Group Theoretic Techniques 

In this appendix, we describe the group theoretic tech- 
niques that we have used to diagonalize the Hamiltoni- 
ans and to study their symmetries. Exact diagonaliza- 
tion [47] is discussed in general in Appendix Al, while 
the symmetries [33, 38] of the irreducible representa- 
tions of -D4 are discussed in Appendix A 2. Finally, in 
Appendix A3, we use symmetry arguments to explain 
why the overall scale of the phase diagrams in the stack 
geometry is larger than the corresponding scale in the 
two-dimensional geometry. Related group theoretic tech- 
niques are discussed in the context of exact diagonaliza- 
tion of the Hubbard Hamiltonian in Refs. [70, 71]. 
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-1 


-1 


s 


Si 


1 
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-1 
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-1 




B 2 
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1 


-1 


-1 


1 


dxy 


E 


2 


-2 











p x and p y 



TABLE III. Character table for the group D4 taken from 
Ref. [47], along with the classification of each representation 
as is discussed in Sec. A 2. 

1. Diagonalization of the Hamiltonian 

As is discussed in Sec. Ill A, the operators n^ and n^ 
commute with the Hamiltonian Eq. (1) and with each 
other, so we diagonalize subspaces with fixed values of 
n^ and separately. The Hamiltonian Eq. (1) has the 
symmetries of a square described by the group -D4 so 
we use the irreducible representations of D4 to further 
simplify the task of diagonalizing each subspace. The 
discussion below uses D4 as an example but is general 
and can be applied to any Unite group. In this section, 
we closely follow Ref. [47]. 

D4 has h = 8 group elements corresponding to the 
symmetries of a square, considered here to be lying in 
the xy plane: the identity E, rotations by 7r around the 
x, y, and z axes C^, C 2y , and C^a, rotations by 7r around 
the lines x = y and x — —y C^ X y and C%xy^ an d coun- 
terclockwise and clockwise rotations by ir/2 around the 
z axis C 4z and Cl}. There are five conjugacy classes 
and thus five irreducible representations for _D 4 . The five 
classes are the identity E, the tt rotation about the z axis 
C2 consisting of C 2z , the two ir/2 rotations 2C4 consist- 
ing of C^z and C^, 1 , the two it rotations about the x and 
y axes 2C2 consisting of C 2x and C 2y , and the two 7r ro- 
tations about the lines x — y and x — —y 2C 2 consisting 
of C 2xy and C 2x y- There are four one-dimensional repre- 
sentations Ai, A 2 , Bi, and B 2 and one two-dimensional 
representation E. The character table for D4 is shown in 
Tab. III. 

The number of times the nth irreducible representation 
appears in the decomposition of a reducible representa- 
tion is [47] 

where the sum is taken over all group elements R, (R) 
is the character of R in the nth irreducible representation, 
and x(R) is the character of R in the reducible represen- 
tation. Let I a) denote a basis function in the occupation 
basis. By writing the eight symmetry operations in the 
occupation basis for a fixed n^ and n^ subspace of the full 
Hilbert space, we calculate the character of each element 
in this occupation representation. Using Eq. (Al), we 
calculate the group structure of each of these subspaces. 
The results are shown in Tab. IV. 
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yields a function 



(n) 



»A(a) 



that transforms as the ith row 



of the nth irreducible representation. For the ith row of 
the nth irreducible representation, a given \a) will only 
be composed of a function belonging to one of the a n 
copies, the A(ct)th copy, of the ith row of the nth irre- 



(n) 
iX(a) 



ducible representation even if a n > 1. Note that 

will be zero if a„ — 0. For each nonzero / , 
responding normalized basis function is, up to a global 
phase, 



fjxL\ ), the cor- 



TABLE IV. Group structure of the single-plaquette Hilbert 
space with D4 symmetry in the presence of n-f- and n_i conser- 
vation. The results are symmetric on interchange of and 
nj.. 



An) 

°iX(a) 



(n) 
i\(a) 



-1/2 



(n) 



iX( a ] 



The coefficients b$ of the decomposition Eq. (A2) are 



As we will see below, the Hamiltonian is block diago- 
nal in the representation basis, so we would like to change 
basis from the occupation basis \a) to the representation 

/ to reduce the dimensions of the matrices 



basis 



dn) 
iX 



refers to a basis 



needing to be diagonalized. Here 

function for the ith row of the nth irreducible represen- 
tation. If the nth irreducible representation is present 
more than once in the decomposition of the occupation 
representation, then there will be a n > 1 orthogonal basis 
functions for the same row of the nth irreducible repre- 
sentation indexed by A. Let there be c irreducible rep- 
resentations and c classes. (For D4, c = 5). \a) can be 

expanded in terms of the 



4 n) 



as 



EEE 6 

Tl=l i=l X = l 



(n) 
iX 



4 n) 



(A2) 



where l n is the dimension of the nth irreducible repre- 
sentation. 

t>ix/ an d are found by applying the projec- 



The 
tion operator 



^ h ^ 

R 



(A3) 



to I a). rw(iZ) is the nth irreducible representation of 
the group element R and Pr is the operator correspond- 
ing to R that acts on functions instead of coordinates 
and satisfies 

P R g{x) = giR-'x). 

To compute Pr in practice, one finds the action of R 
on the spatial indicies of the creation operators used to 
define the second-quantized wave function. 

(n) 

The projector projects into the ith row of the nth 
irreducible representation, so 



iX(a) 



Pl n) la) 



(A4) 



6 («) _ 
u iX ~ 



if a n = 

p(") I f («) 



1/2 



(flxU\mU)) 5 A)A(Q) ifa n ^0 



To find all of the basis functions for the representation 
basis, we compute Eq. (A4) for every basis function \a) 
in the occupation basis for all P^ [72]. By applying 
the same projection operator to different functions |cv), 
different functions for the same row of the same repre- 
sentation will be generated for all A = 1, ... , a n . Once 



all of the 



An) 
} iX 



have been found, the Hamiltonian can 



be transformed into the representation basis by 

-^representation — & -^occupation S : 

where S is the transformation matrix given by S a ^ = 
(al^,-^). Since Pr commutes with H for all i?, 



} iX(a) 



\H\ 



Am) 



p( n ) frp( m ) 

a 33 



P^j OC 8 n , m Si 



Thus, matrix elements between different irreducible rep- 
resentations or different rows within the same irreducible 
representation vanish. Therefore, this transformation 
into the representation basis can be made separately for 
each row of each irreducible representation, and each 
of the resulting Hamiltonians can be diagonalized sepa- 
rately. This greatly reduces the dimension of the Hamil- 
tonians that need to be diagonalized since the dimension 
of the Hamiltonian for the ith row of the nth irreducible 
representation is a n . States found in the representation 
basis can be transformed back into the occupation basis 
by 

(<#>=$3S a ,*<^>. 

4> 

In order to compute the projector Eq. (A3), it is 
necessary to know the representation T^ n \R) of each 
group element R in each representation n. For the one- 
dimensional representations, the representations T^ n \R) 
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are just the characters x^ n \R) listed in Tab. III. In the 

x 



basis 



the representations of the group elements in 



the E representation are 



1 
1 



T {E) {C 2x ) = 



-1 o 
-1 



-1 

1 

1 
-1 

1 

1 



r [E) (c 2xy ) = 



o 1 

-1 

-1 
1 

-1 
-1 



(A5) 



2. Classification of the Eigenstates 



In this section, we classify the symmetries of the five 
representations of D4 and give examples of states trans- 
forming as each of these representations. Since the 
A and B representations are one-dimensional, we can 
study their symmetries directly from the character table 
Tab. III. The A states are symmetric under ir/2 rotations 
2C4, while the B states are antisymmetric under these 
rotations. For this reason, the A states are classified as 
s-wave and the B states are classified as d-wave. The Ai 
representation is symmetric under all five classes, while 
the A 2 representation is antisymmetric under 2C 2 and 
2C 2 . The Bi representation is symmetric under 2C' 2 and 
is antisymmetric under 2C 2 , while the B 2 representation 
is antisymmetric under 2C 2 and is symmetric under 2C 2 . 
Thus, B\ is classified as d x 2_ y 2 and B 2 is classified as d xy . 

The E representation is two-dimensional, so we need 
to consider the representations of the elements given by 
Eq. (A5). Both rows of the representation are antisym- 
metric under 7r rotations about the z axis. For this rea- 
son, E is classified as p-wave. The first row transforms 
into the second row and the second row transforms into 
negative the first row under a positive rotation by tt/2 
about the z axis. The first row is antisymmetric and the 
second row is symmetric under a n rotation about the x 
axis. The first row transforms into the second row and 
the second row transforms into the first row under a tt 
rotation about the line y = x. Similar considerations for 
the other rotations lead to the classification of the first 
row of the E representation as p x and the second row as 
Pv 

Simple examples of states in the A\, B 2 , and E repre- 
sentations come from n+ = 1 and — 0. An example of 
an Ai state is 



1 



We use the notation in which, for example, 
It I 0) = 4 t c^ |0> , 

and the numbering of sites within a plaquette is given in 
Fig. 1. This state is clearly invariant under all five of the 
D4 classes. An example of a B 2 state is 



\B*) = \ (It 000) 



|ot 00) + |ooto) - |000 t». 



This state is invariant under C 2 and 2C 2 but changes 
sign under 2C4 and 2C' 4 . An example of two E states is 

\E X ) = \ (It 000) -|0t 00) - [00 t 0) + |000 t)) , 

\E V ) = \ (It 000) + |o 1 00) - [00 to) - 1000 t)) • 

These states are antisymmetric under 7r rotations about 
z. Under a positive rotation by tt/2 about z, \E X ) — > \E y ) 



and \E„ 



— \E X ). Under a tt rotation about the x axis, 



\E X ) — > \E y ) and \E y ) — ► — \E X ). Under a n rotation 
about the line y — x, \E X ) — > \E y ) and \E y ) — » \E X ). 
Thus \E X ) transforms as p x and \E y ) transforms as p y . 

A simple example of an A 2 state comes from n-t- = 2 
and = where 



\M) = 2 (Itt 00) + |0 tt 0) 



1 00 tt) - It 00 1)) 



2 Y c rf c 2t 



L 3t 4f 



°4t L lt 



(It 000) + |0 1 00) + 1 00 1 0) + |ooo t)) ■ 



This state is invariant under C 2 and 2C4 since this cor- 
responds to cyclically permuting the indices. Under 2C' 2 
and 2C 2 , \A 2 ) — > — \A 2 ) since this corresponds to swap- 
ping indices. 

A simple example of a B\ state is just the state from 
— 4 and = 

\Bi) = itttt) = 4 t 4t4r4t i°) ■ 

As can be seen by appropriately switching the indicies on 
the creation operators, this state has the B\ symmetries 
discussed above. 



3. Scale Increase in the Phase Diagrams for the 
Stack Geometry 

In this section, we use symmetry arguments to explain 
the increase in the overall scale of phase diagrams for the 
stack geometry (Figs. 10, 11, and 12) relative to the cor- 
responding diagrams for the two-dimensional geometry 
(Figs. 3, 7, and 8). The main reason for this increase 
is that the J' ± , J' z , and W transition matrix elements 
between the low- and high-energy subspaces decrease as 
J" I J' and J'" I J' are varied from to 1. As a result, all 
phase transitions occur at smaller values of t 1 . To gain 
intuition for why these matrix elements decrease as J" 
and J'" approach J', we study the t-J± Hamiltonian for 
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the stack geometry in the regime where J' ± = J'J_ = J"' 
(obtained as a' /a — > co). A similar argument holds for 
the other Hamiltonians considered. 

Consider two neighboring plaquettes in a stack in the 
1 2 , 2} 1 2, 2} state coupled by the J' ± perturbing Hamilto- 
nian H 1X with J' ± = J'l = Jf . From Tab. II, both 
plaquettes have B\ symmetry. Since every vertex in one 
of the plaquettes is equally coupled to every vertex in the 
other plaquette, H±± is symmetric under arbitrary per- 
mutations of the vertices within one of the plaquettes. 
As a result, the matrix element 

(l,3|<3,l|ff 1± |2,2)|2,2) 

will vanish unless both |1, 3) and |3, 1) have B\ symmetry; 
otherwise, an appropriate permutation can be used to 
show that this matrix element is equal to negative itself. 
This drastically reduces the magnitude of /j 4 ' 4 ^ relative 
to its value in the J" = J"' = case. A similar argument 



shows that /j 4 ' 2 ^ 



and /j 2,2 ** vanish for J' ± 



77/ _ Till 

(since the |1,3) and |3, 1) manifolds have no A\ states. 
See Tab. IV). 



Appendix B: Basis Vectors in the Representation 
Basis 

In this appendix, we list the basis vectors for the irre- 
ducible representations for 7D4. The results are symmet- 
ric on interchange of up spins with down spins. We use 
the notation |n^,n^(r)) to refer to the basis vector for 
a plaquette with rif up spins and down spins in the 
irreducible representation V. If there are a n > 1 copies 
of an irreducible representation, we denote them by a su- 
perscript. We use subscripts x and y to denote the first 
and second rows of the E representation, respectively. 
Thus |l, l(Ey)} refers to the basis vector for |1, 1) in the 
first copy of the second row of the E representation and 
|2, l(Bf)') refers to the basis vector for |2, 1) in the second 
copy of the B\ representation. 
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Appendix C: Explicit Expressions for Important 
Ground States and Energies 

In this appendix, we give explicit expressions for the 
ground states and energies for those |2) and |4) states 
responsible for d-wave superfluidity in the t- J± Hamilto- 
nian, t-J±-J z Hamiltonian, and the t-J-W Hamiltonian 
with Jj_ = J z = J. We also discuss intuition for the 
symmetries of some of these states. 



1. t-J± Hamiltonian 

Here we give the |2) and |4) ground states of the t- 
J± Hamiltonian responsible for <i-wave superfluidity (see 
Tab. II and Fig. 2). The s-wave symmetric |1, 1) ground 
state for J±/t > —1.22 is 



|l ! l(A 1 ))oc6|l ! l(A 1 )) + |l,l(A?)) 
with energy 



(CI) 



£ 9 (M) 



1 
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In- N\ '2 j ( ^ +2048 



(4 + 72) 



J± 
t 



(C2) 



where 
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1024 +(2V2- 1) 



J \ 



The d-wave symmetric |2, 2) ground state for Jj_/t > 



is 



./fie: _ 1 

12,2(50) « ^— - \2,2(Bl)) - \2,2(Bf)) (C3) 



with energy 



E g (2,2) 



65 J ± 



4V2 t 



We can understand the symmetries of the |2) ground 
states as follows. In the limit of vanishing t, the ground 
state will be an eigenstate of the J± interaction. The two 
eigenstatcs of the Jj_ interaction for an up and a down 
molecule on sites r and r' axe 



\±) 



^2 ( C U°r'l ± C ti C t't 



(C4) 



^ (s+s-, + s-s+) \±) rr , = ±^ |±)„ 



(Note that 



sl r , |0) where s T , is given by 



Eq. (4)). Thus for positive Jj_, the singlet state 
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is the ground state of the J± interaction. Note that it 
is disadvantageous to have identical spins since, in that 
case, the J± interaction vanishes. For two molecules with 
small tunneling amplitude t, the |2) ground state is thus 
a superposition of singlets on the four nearest-neighbor 
bonds. Diagonalizing the effective Hamiltonian of this 
system shows that the ground state is the following sym- 
metric superposition: 



/23 



'34 



This state has s-wave symmetry and is proportional to 
the basis vector |1, l(Ai)), consistent with Eq. (CI). For 
negative J±, the |+) state is the ground state of the J± 
interaction on two sites, and the two resulting single- 
plaquette ground states have p-wave symmetry. 

The symmetries of the |4) ground states can be under- 
stood by considering only nearest-neighbor interactions. 
In this case, for positive Jj_, the ground state is 

|2,2(^)>-|2,2(B 1 1 )>, 

which has d x 2_ y 2 symmetry. For negative J]_, the ground 
state is 

|2,2(^)) + |2,2( J B 1 1 )), 

which also has d x 2_ y 2 symmetry. When the effects of 
next-nearest neighbors are considered, the coefficients 
in the superposition of |2, 2{B\)) and |2,2(_B^)) are 
changed. 



2. t-J±-J z Hamiltonian 

Here we give the |2) and |4) ground states of the t-J±- 
J z Hamiltonian responsible for d-wave superfluidity (see 
Fig. 4). The s-wave symmetric |1, 1) ground state is 

|l,l(A 1 ))(xa|l,lK)) + |l,l(A?)) (C5) 

with energy 
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16 



where 
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The ci- wave symmetric |2, 2) ground state is 

|2,2(Bi))a6|2,2(Si)) + |2,2(B?)) (C6) 



with energy 



E g (2,2 




2^-16^ 
t t 



where 




„i^) 2 +8 ( 2 ^-i 



1/2 



+ t + ( 8 ^- 4 )t 



3. t-J-W Hamiltonian with J± = J z = J 

Here we give the |2) and |4) ground states of the t-J- 
W Hamiltonian with J± = J z = J responsible for e?-wave 
superfluidity (see Figs. 5 and 6). The s-wave symmetric 
|1, 1) ground state is 

|l,l(A 1 ))(xa|l,lK)) + |l,l(A?)) (C7) 

with energy 



E g (l,l) = ~ 



18 9 - 4V2 



2048+ 



(4 + ^2) 



(C8) 



where 
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The d-wave symmetric |2, 2) ground state is 

12,2(5!)} tx b \2,2(Bl))+\2,2(B*)) (C9) 
with energy 

1 J " 



E g (2,2) 



where 



b = - 



V2-16) - V402 - 96V2 



>2- 



V201 - 48\/2 J 



t 



We can understand the symmetries of the ground 
states of this Hamiltonian in the different sectors as fol- 
lows. For large positive W, when the ground states are 
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|0, 2) and |0, 4), the |4) Hilbert space is one-dimensional 
with state 

10,4(5!)) = |||||), 

so the 1 4) state automatically has d x 2_ y 2 symmetry. In 
the limit of infinite W, the |2) ground state is a superpo- 
sition of states with two down spins next to each other. 
In this limit, the states in the A 2 and B 2 representations 
and the two states 

|0, 2{E X )) ex |0 || 0) + || 00 1> + O (J^j , 
|0, 2{E y )) <x |00 ||) + HI 00) + O (j^j 

in the E representation are the ground states. For fi- 
nite but large W, the 0{t/W) terms, proportional to 
|0 | |) and || | 0), make the p-wave symmetric states 
the ground states. 

As W is decreased, the |2) ground state is s-wave 
|1, 1) since positive J± and J z make it favorable to have 
molecules of opposite spins. For intermediate values of 
these interaction strengths, the |4) ground state becomes 



1 1,3). With four molecules on a plaquette, there is no 
tunneling to consider in the energies of the states. Fur- 
thermore, since and are fixed, as determined by J z , 
Jj_, and W in this sector, we only have to consider the 
effects of J± to explain the behavior of |4) in this sector. 
Consider two molecules of opposite spins on sites r and 
r' . The two eigenstates of the J± interaction are |±) . 
defined in Eq. (C4). Since J± > 0, the singlet |— ) rr , is 
the ground state of the J± interaction. Therefore, we 
may expect the |1, 3) ground state to be a superposition 
of singlets on nearest-neighbor bonds with relative phases 
chosen for constructive interference: 




Here 

t rr ' — c\^c\,^ 

creates the m = — 1 triplet state on sites r and r', while 
s rr i , defined in Eq. (4) , creates the singlet state on sites 
r and r'. This state is indeed the |1, 3) ground state and 
has the s-wave A 2 symmetry. 
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